Persistence to high temperatures of interlayer coherence in an organic superconductor 
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The interlayer magnetoresistance p zz of the organic metal k-(BEDT-TTF)2Cu(NCS)2 is studied 
in fields of up to 45 T and at temperatures T from 0.5 K to 30 K. The peak in p zz seen in in- 
plane fields, a definitive signature of interlayer coherence, remains to Ts exceeding the Anderson 
criterion for incoherent transport by a factor ~ 30. Angle-dependent magnetoresistance oscillations 
are modeled using an approach based on field-induced quasiparticle paths on a 3D Fermi surface, 
to yield the T dependence of the scattering rate r _1 . The results suggest that r _1 does not vary 
strongly over the Fermi surface, and that it has a T 2 dependence due to electron-electron scattering. 
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The past two decades have seen a tremendous blos- 
soming of interest in compounds that possess quasi-two- 
dimensional (Q2D) electronic bandstructure; examples 
include crystalline organic metals 0, 0, S 3 1 cuprates Q 
and layered ruthenates @- These materials may be de- 
scribed by a tight-binding Hamiltonian in which the ratio 
of the interlayer transfer integral t\ to the average in- 
tralayer transfer integral t\\ is <C 1 [IIlJlJlI The ques- 
tion arises as to whether the interlayer charge transfer is 
coherent or incoherent in these materials, i.e. whether 
or not the Fermi surface (FS) is three dimensional (3D), 
extending in the interlayer direction. Various criteria for 
interlayer incoherence have been proposed, including 



k B T > t±, 



(1) 



where T is the temperature. In this picture, thermal 
fluctuations "wipe out" details of the interlayer peri- 
odicity 7]. The consequent interlayer incoherence is 
used as a justification for a number of theories which 
are thought to be pivotal in the understanding of Q2D 
materials (see e.g. 0, H, El It is therefore im- 

portant to test assertions such as Eq. ^ To this end, 
we have measured the magnetic-field-orientation depen- 
dence of the resistance of the organic metal re-(BEDT- 
TTF) 2 Cu(NCS) 2 @ using fields of up to 45 T and Ts of 
up to 30 K (Fig-QJ. This material was chosen because its 
FS is well knownpllllj (see Fig.^c)), previous low-T ex- 
periments have shown that the interlayer transfer integral 
is f a = 0.065 ± 0.007 meV [HEll and modest Ts allow 
the inequality in Eq. [2 to be exceeded by orders of mag- 
nitude (i a /&B w 0.5 K). Our data show that interlayer 
coherence survives to at least k^T >30i a (c.f. Eq.QJ. In 
addition, the results in this paper determine the low-T 
scattering rate r _1 for the quasiparticles. Recent works 
stress that an understanding of the effects of the scatter- 
ing mechanism in organic metals gives vital information 
about the mechanism for superconductivity 0, Q] . 
Electrical contacts were made to single crystals (~ 



0.7 x 0.5 x 0.1 mm 3 ) of k-(BEDT-TTF) 2 Cu(NCS) 2 us- 
ing 12.5 /itm Pt wire bonded by graphite paste; contact 
resistances were ~ 10 Q. Current and voltage termi- 
nals were arranged such that the measured resistance 
R zz is proportional to the interlayer component of the 
resistivity tensor, p zz The measurements in Figs.Q] 
and [21 use a well-characterized crystal also employed in 
Rets^! El El- Additional experiments (e.g. Fig. 0) 
used d8-K-(BEDT-TTF) 2 Cu(NCS) 2 crystals; here, "d8" 
indicates that the terminal hydrogens of BEDT-TTF 
have been replaced by deuterium. These crystals were 
studied in Ref. £ 11], where it was found that deuteration 
reduces t a to 0.045 ± 0.005 meV. Samples were mounted 
on a ceramic holder attached to a cryogenic goniometer 
allowing continuous rotation in static magnetic fields B. 
Crystal orientations are labeled by the angles 0, <fi. Here, 
9 is the angle between B and the normal to the highly- 
conducting planes and <fi defines the plane of rotation; 
4> = is a plane of rotation containing the fcb direc- 
tion and the normal to the highly-conducting planes . 
Sample Ts were stabilized using a calibrated Cernox sen- 
sor mounted on the sample holder and a heater driven 
by a feedback circuit. The Cernox resistance was mon- 
itored during isothermal field sweeps in order to deduce 
the correction to the apparent T due to magnetoresis- 
tance; results were in good agreement with the study of 
Ref. [l6j ■ Static fields were provided by the 45 T Hybrid 
magnet at NHMFL Tallahassee and 33 T Bitter magnets 
at Tallahassee and HFML Nijmegcn. 

Typical R zz data are shown in Figs, ^a) (4> — 160°) 
andQIb) (0 = 90°) for Ts in the range 5.3 K to 29.3 K. 
In the former rotation plane, features in p 7Z are known 
to be due to orbits on the Q1D FS sections whereas 
in the latter they can be confidently attributed to Q2D 
FS orbits A series of angle-dependent magnetore- 

sistance oscillations ( AMROs) , periodic in tan 9 |9j, [llj is 
observed at both (j>; those due to the Q2D FS sections are 
frequently known as Yamaji oscillations [9j . In addition, 
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FIG. 1: Interlaycr resistance R zz (oa p zz ) of a k-(BEDT- 
TTF)2Cu(NCS)2 crystal versus tilt angle 9 for various con- 
stant T; B = 45 T. (a) Data for <f> = 160°, a plane of rotation 
at which p zz is determined by phenomena on the Q1D FS 
sections. In order of increasing R zz at 9 = 35°, the curves are 
for T = 5.3, 6.5, 7.6, 8.6, 9.6, 10.6, 12.1, 13.1, 14.6, 17.1, 19.6 
and 29.3 K respectively. In addition to AMROs, DKCOs can 
be seen as twin peaks on either side of 9 = 90° in the lower 
T data, (b) Similar data for <f> = 90°; here p zz features are 
associated with the Q2D FS section. In order of increasing 
R zz at 9 = -70°, the curves are for T = 5.3, 7.6, 8.6, 9.6, 
10.6, 12.1, 13.1, 14.6, 17.1, 19.6, 24.5 and 29.3 K respectively. 

Danner-Kang Chaikin oscillations (DKCOs) [Til Il7| are 
observed as peaks either side of 9 = ±90° when <\> — 160°. 
As T increases, the AMROs decrease in intensity, with 
the higher-index oscillations (i.e. those closer to 9 = 
90° [lj) disappearing first. After most AMROs disap- 
pear, the slowly- varying background magnetoresistance 
begins to increase more markedly. 

FiguresEJa) and (b) show expansions of data in Fig.Q] 
close to 9 = 90° at selected Ts. For ease of compari- 
son, data are plotted as Ap zz / p ZZ BG, the fractional de- 
viation in p zz from the background magnetoresistance 
PzzBG determined by fitting a smoothly- varying curve 
through the data on either side of 9 — 90° Q|. In 
such measurements, interlayer coherence is detected us- 
ing a phenomenon known as the "coherence peak" or 
"SQUIT (Suppression of QUasiparticle Interlayer Trans- 
port) peak" p], Q) lJ M- LUj > a maximum in Ap zz / p ZZ BG 
observed when B lies exactly in the intralayer plane (i.e. 
at 9 = 90°). This occurs because of the efficient inter- 
layer velocity averaging caused by closed orbits on the 
side of the FS (Fig|2Ic-e)); these can exist if, and only 
*/QjBi the interlayer transport is coherent, i.e. the FS is 
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FIG. 2: (a) The 45 T magnetoresistance close to 9 = 90° 
at (f> = 160°, plotted as Ap zz /p ZZ BG, the fractional change 
in p zz from the more slowly-varying background. Data for 
T = 5.3 (highest), 7.6, 9.6 and 13.1 K (lowest) are shown, 
offset for visibility, (b) Similar data for <f> — 90°; the traces 
are for T = 5.3 (highest), 7.6, 8.6, 9.6, 10.6, 12.1, 13.1 and 
14.6 K (lowest). Each trace has been offset for clarity, (c) 3D 
representation of the FS of k- (BED T-TTF) 2 Cu (N C S ) 2 ( aft er 
Ref. the finite t a gives the corrugations (shown greatly 

exaggerated) on the sides of the FS. Q1D and Q2D FS sections 
are shown in red and blue respectively, (d) Consequent field- 
induced closed orbits on the side of the Q1D FS section when 
9 = 90° and = 0. (e) Similar closed orbits on the Q2D FS 
section when 9 = 90° and (j> — 90°. Orbits such as those in 
(d) and (e) give rise to the SQUIT peak in p zz . 

a 3D entity extending in the interlayer direction. In both 
Figs. EJa) and (b), the SQUIT peak is visible close to 
9 = 90°. In spite of the small size of t a , this demonstra- 
tion of interlayer coherence continues to be observable up 
to at least 14.6 K, exceeding the criterion in Eq. ^by a 
factor ~ 30 p|. 

The AMROs and DKCOs in Fig. 2] are very similar to 
those observed at T = 0.5 K [Tj, the only difference be- 
ing that they decrease in amplitude as T increases; i.e. 
the magnetoresistance features do not change in form or 
angular position as T grows, they merely fade gradually. 
This suggests that the same mechanism is responsible 
for the form of the magnetoresistance at all T examined. 
Moreover, the SQUIT peak (Fig. |2Jl demonstrates that 
the FS of k-(BEDT-TTF) 2 Cu(NCS) 2 remains 3D up to 
at least T « 15 K 18]. We therefore choose to simu- 
late the data using the numerical method based on field- 
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FIG. 3: Comparison of experimental p zz data (a) and the nu- 
merical simulation (b) method described in Ref. [Tl|. T = 
1.5 K and <j> = 15°; for the lowest to the highest traces, 
B is 32, 34, 36, 38, 40, 42 and 44 T (no offset is applied), 
(c) AMRO data for <f> = 150°, T = 1.5 K and fields of 20 
(lowest), 24 and 28 T (highest). The arrow indicates the am- 
plitude AR ZZ of a particular AMRO feature; < R zz > is the 
resistance at the arrow's midpoint, (d) Experimental AMRO 
AR ZZ / < R zz > for <f) = 150° plotted as a function of the 
orbit angular frequency ui times scattering time r. Data for 
several AMROs, denoted by their Yamaji indices Q i are 
shown (i = 2 (square), i — 3 (dot), i = 4 (triangle) and i — 5 
(inverted triangle)) for several Ts in the range 1.7 K to 5.5 K. 

induced quasiparticle trajectories across a 3D FS that 
was used to model T = 0.5 K AMRO results successfully 
in Ref. ^lj. This procedure employs Chambers' Equa- 
tion and the experimentally-deduced FS (Fig. |2{c) 11]); 
given B and the angles 9, cj>, the only input parameter re- 
quired to calculate p zz is r. The model has the restriction 
that magnetic breakdown between the Q1D and Q2D FS 
sections [2(j is not included; recently it was shown that 
breakdown in k-(BEDT-TTF) 2 Cu(NCS) 2 can produce a 
series of AMROs [2lJ. However, the breakdown ampli- 
tude falls off rapidly with increasing 9 [2(j , so that these 
effects are small for B rj 45 T when 9 > 70° In the 
following, we neglect data at lower 9 for this reason. 

A typical comparison of experiment and simulation is 
shown in Figs. Ola) and (b) (<fi — 15°). For each <f>, 
it is possible to map the normalized AMRO amplitude 
&Rzz/ < Rzz > (see Fig.|^c)) from the simulation onto 
the experimental value by adjusting r. Moreover, for 
each T, it was found that the same r, to within experi- 
mental errors, fitted all the AMRO features. Given this 
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FIG. 4: Scattering time r deduced from AMROs at <f> « 
160° (a) and (f> « 90° (b) versus T. Consistent values of r 
are deduced from several AMRO features at fields of 30 T 
and 45 T (insets show field and 9 for each feature). The data 
are fitted to r = (r^ + AT 2 )' 1 (curves) with to = 3.3±0.1 ps 
and A = 0.0062 ± 0.0003 ps^K -2 (a) and To = 3.5 ± 0.1 ps 
and A = 0.0065±0.0003 ps" 1 !^ 2 (b). These scattering rates 
are identical to within experimental errors, even though the 
AMROs in (a) and (b) are produced by different FS sections. 

result, it is useful to find a single parameter to describe 
all AMRO amplitudes for a particular (j>. A good can- 
didate is lot, where uj is the angular frequency at which 
the orbit responsible for the AMRO feature is traversed. 
In the case of AMROs due to the Q2D FS sections, 
lu = {eB/m* a ) cos 6*, the cyclotron frequency 0. Here, 
m* is the 9 = effective mass of the a Q2D pocket 
of the Fermi surface 0, |2£j. For orbits on the Q1D 
Fermi surface, the relevant period is the time to cross 
the Brillouin zone [23|. As long as one avoids the region 
close to 9 = 90° where the SQUIT and DKCOs occur 
(the latter involving orbits that "snake" across contours 
on the FS 0,0]), this frequency can be expressed as 
lu = {2eB/{m*p - m* )) cos 9, where is the 9 = effec- 
tive mass of the " /?" breakdown orbit [2(j . 

Using these us, the Ai? zz / < R zz > values of all of 
the AMRO features at a particular <f> collapse onto a sin- 
gle curve as a function of ut; an example is shown in 
Fig-Eld) for Yamaji oscillations due to the Q2D FS. Once 
this correspondence is known for each cj), t can be read- 
ily found from AMRO data. The r values thus obtained 
follow a T-dependence of the form r _1 = t^ 1 + AT 2 
rather closely. Moreover, to within experimental errors, 
the scattering rates have the same tq and values of A, 
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irrespective of the FS orbits involved. This is illus- 
trated in Fig. 01 which shows r deduced from AMROs 
at (f> = 160° ((a)- due to Q1D FS sections) and cj> = 90° 
((b)- due to Q2D FS sections). In spite of the fact that 
the quasiparticle orbits involved are very different, and 
involve distinct FS sections, the T dependence of r for 
both is virtually identical. This suggests that mecha- 
nisms for superconductivity in organic metals that in- 
voke a large variation in scattering rate over the FS {e.g. 
"FLEX" methods 0) are inappropriate for k-(BEDT- 
TTF) 2 Cu(NCS) 2 . For all FS orbits studied, the inferred 
T = scattering time (r = 3.4 ± 0.2 ps- see Fig. gj) is 
close to values measured by other means [T5I- 

The T 2 dependence of r _1 suggests electron-electron 
scattering [24| . A T 2 scattering rate was inferred from 
B = resistivity measurements |25j . However, prob- 
lems in deconvolving the in-plane resistivity component 
P|| from p zz in experimental data [fit Il5|. and the influ- 
ence of the broad superconducting transition on the In- 
dependence of the measured resistivity [jj mean that this 
attribution can not be considered conclusive. By con- 
trast, the T-dependent AMROs provide an unambiguous 
gauge of the scattering rate of normal-state quasiparti- 
cles, allowing the mechanism to be definitively identified. 

When k-gT > 4i a , the interlayer contribution to the 
bandwidth, the FS shown in Fig. [2fc) is "blurred" on a 
wavevector scale k^T/hvp (where vp is the Fermi veloc- 
ity) exceeding the amplitude of interlayer corrugation. 
Nevertheless, a Fermi-liquid picture predicts that the 
semiclassical dynamics of each electron are still governed 
by an equation of motion which contains and reflects the 
detailed nature of the interlayer dispersion. Quasiparticle 
states within ~ k^T of the chemical potential continue 
to contribute to the AMROs and SQUIT, as long as the 
underlying dispersion is not affected by a variation in 



electron energy on the scale of k^T; this will hold when 
k B T < t|| ~ 10 - 100 meV |ll|. Then, the only effect of 
raising T is from T-dependent contributions to the scat- 
tering rate. This is exactly what we find in our AMRO 
data, and demonstrates that a coherent 3D FS picture is 
valid for this compound, at least at T<0.1i|| (a limit our 
experiments do not exceed), at which point a crossover 
to an incoherent regime might be expected j2(|. 

In summary, the orientation dependence of the 
interlayer magnetoresistance p zz of k-(BEDT- 
TTF) 2 Cu(NCS) 2 has been studied in fields of up 
to 45 T as a function of temperature T. The "SQUIT" 
or coherence peak in p zz seen in exactly in-plane fields, a 
definitive signature of a 3D Fermi surface (i.e. interlayer 
coherence), persists to temperatures exceeding the 
proposed Anderson criterion for incoherent transport 
(Eq. 1 7j) by a factor ~ 30. Features in the magne- 
toresistance have been successfully modelled using an 
approach based on field-induced quasiparticle paths on 
a 3D Fermi surface, to yield the temperature depen- 
dence of the quasiparticle scattering rate. The results 
suggest that the average scattering rate does not vary 
strongly over the Fermi surface, and that it follows a T 2 
dependence attributable to electron-electron scattering. 
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